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Abstract 

The type-I region of phase transitions at finite temperature of the U (l)-Higgs theory 
in 3+1 dimensions is investigated in detail using a Wilsonian renormalisation group. 
We consider in particular the quantitative effects induced through the cross-over of 
the scale-dependent Abelian charge from the Gaussian to a non-trivial Abelian fixed 
point. As a result, the strength of the first-order phase transition is weakened. Ana- 
lytical solutions to approximate flow equations are obtained, and all characteristics 
of the phase transition are discussed and compared to the results obtained from per- 
turbation theory. In addition, we present a detailed quantitative study regarding the 
dependence of the physical observables on the coarse-graining scheme. This results 
in error-bars for the regularisation scheme (RS) dependence. We find quantitative 
evidence for an intimate link between the RS dependence and truncations of flow 
equations. 
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I. INTRODUCTION 



The phase transition of the {/(l)-Higgs theory in 3+1 dimensions at finite temperature pro- 
vides an important model for cosmological phase transitions. In the high temperature limit, 
it reduces to the purely 3d Abelian Higgs model describing the superconducting phase tran- 
sition |l| , or certain nematic to smectic-A phase transitions in liquid crystals . The phase 
transition in this model is governed by the infra-red region of its spectrum of fluctuations. 
The nature of the phase transition depends primarily on the ratio mn/mA between the scalar 
and the gauge field mass. For superconductors, these mass scales correspond to the inverse 
correlation length and the inverse London penetration depth, respectively. For small values 
of the Higgs field mass the phase transition is strongly enough first order to cut-off long 
range fluctuations. This corresponds to the good type-I region for standard superconduc- 
tors, mn/mA < 1- On the other hand, the type-II region corresponds to mn/m^ > 1. Here, 
it is expected that the phase transition changes from first to second order. 

A proper treatment of the long range fluctuations is decisive for an understanding of the 
U (l)-Higgs phase transition as they change the effective interactions between the fields. The 
'microscopic' physics in the ultra-violet is characterised by the couplings at the short dis- 
tance length scale 1/A (or ~ 1/T, with T the temperature). In turn, the physics close to 
a first-order phase transition depends typically on the (small) photon mass rriA, and thus 
requires the knowledge of the couplings at scales <C T.[] 

A field-theoretical approach which in principle is able to deal with the effects of long range 
fluctuations and which describes the related scaling of the couplings is given by the Wilsonian 
renormalisation group . This procedure is based on integrating-out infinitesimal momen- 
tum shells about some 'coarse-graining' scale k within a (Euclidean) path- integral formula- 
tion. The infra-red effective theory obtains upon integrating the resulting flow w.r.t. k — > 0. 
This way, the characteristic scaling behaviour (or 'running') of the couplings as functions 
of k, and in particular the running of the Abelian charge e(k), is taken into account. A 
Wilsonian approach thus improves on perturbative resummations in that the perturbative 
expansion parameter e 2 T/rriA becomes now scale dependent, e 2 (k)T/rriA{k). While the for- 
mer diverges close to a second-order phase transition where the photon mass vanishes, the 
latter remains finite in the infra-red even for rriA{k) —>■ due to the non-trivial scaling of the 
Abelian charge. The crucial role of running couplings in finite temperature phase transitions 
has been discussed in pure scalar theories [^,|[. 

In the present article we employ the Wilsonian renormalisation group to the type-I regime 
of the [/(l)-Higgs phase transition. Our main contributions are two-fold. First, we take 
into account the non-trivial scaling of the Abelian charge e 2 (k), characterized by an effec- 
tive Abelian fixed point which is kept as a free parameter. The infra-red effects lead the 



1 When speaking of 'scales' we have always 'mass scales' or 'momentum scales' ~ k in mind. The 
corresponding 'length scales' are given as ~ k^ 1 . 
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Abelian gauge coupling to cross-over^ from its slow logarithmic running in the ultra-violet 
(effectively 4c?) to a strong linear running in the infra-red (effectively 3d). The characteris- 
tic scale for this cross-over depends on the precise infra-red (IR) behaviour of the Abelian 
charge, and is decisive for both the strength of the transition and the properties of the 
phase diagram. This is currently the least well understood part of the problem. Equally 
important is to retain the full field dependence of the effective potential (no polynomial 
approximation), for which an analytical expression is given in the sharp cut-off case. We 
obtain all thermodynamical quantities related to the first-order phase transition and study 
their dependence on the cross-over behaviour. Second, we present a detailed quantitative 
analysis of the 'coarse-graining' dependence of our results. This is an important consistency 
check for the method and the approximations involved. We give quantitative evidence for 
an intimate link between a truncation of the effective action, and the dependence on the 
coarse graining scheme, which can simply be displayed as additional 'error-bars' due to the 
scheme dependence. 

The 3d [/(l)-Higgs phase transition has been studied previously using flow equations PHT3 



and within perturbation theory flT4}4T6||. Recent results from lattice simulations for both 



type-I and type-II regions have been reported as well [17,18]. In M, the RG flow has been 



studied for the type-II regime within a local polynomial approximation for the effective po- 
tential about the asymmetric vacuum up to order 8 in order to establish the phase diagram, 
the relevant fixed points and the related critical indices. The polynomial approximation is 
expected to give reliable results for the scaling solution close to a second-order fixed point. 
The type-I regime has been discussed for the full potential, using a matching argument for 
the running of the Abelian charge. In ||10|| , the large -N limit and its extrapolation down to 
N = 1 has been considered as well. It was pointed out that the local polynomial approxima- 
tion becomes questionable close to a first-order phase transition or a tri-critical fixed point 
at about N ~ 4 or smaller. This was later confirmed by Tetradis |12|] , who in addition aban- 
doned the local polynomial approximation. The present study, aiming particularly at the 



type-I region of the phase diagram, improves on P,|T0|| in that the full field dependence of the 
effective potential will be taken into account. A quantitative description of thermodynam- 
ical observables at the phase transition requires a good accuracy for the effective potential 
in the first place. Our study also goes beyond the work of [12[] in three important aspects. 
We study the dependence of physical observables on the value of the effective Abelian fixed 
point. In addition, explicit analytical solutions to approximate flow equations are given, as 
well clS db discussion of the scheme dependence. 

This article is organised as follows. We introduce the Wilsonian flow equations and the 
particular Ansatz used for the £7(1)-Higgs theory. The flows for the Abelian charge and 
the free energy are explained, as well as the further approximations involved (Sect. |T|). We 
then proceed with the thermal initial conditions as obtained from perturbative dimensional 
reduction (Sect. pip and a discussion of the phase diagram and the critical line (Sect. [TV| ). 
This is followed by a computation of all relevant thermodynamical quantities at the first- 



2 This cross-over is not to be confused with the qualitatively different 'cross-over' observed in the 
type-II regime of 3+1 dimensional SU (2) +Higgs theory. 
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order phase transition as functions of the effective Abelian fixed point, a computation of the 
corresponding characteristic scales, and a discussion of the approximations made (Sect. 0). 
A quantitative study of the scheme dependence on the main characteristics of the phase 
transition is given (Sect. [VI]), followed by a summary and an outlook (Sect. [VI H ). Three 
Appendices contain some more technical aspects of our analysis. 



II. FLOW EQUATIONS 
A. Wilsonian flows 

Wilsonian flow equations are based on the idea of a successive integrating-out of momen- 
tum modes of quantum fields within a path-integral formulation of quantum field theory 
Hf§]. This procedure, in turn, can also be interpreted as the step- by-step averaging of the 
corresponding fields over larger and larger volumes, hence the notion of 'coarse- graining'. 
The modern way of implementing a coarse-graining within a path-integral formalism goes 
by adding a suitable regulator term ~ / 0i?fc(g)0, quadratic in the fields, to the action [Q. 
This additional term introduces a new scale parameter k, the 'coarse graining' scale. A 
Wilsonian flow equation describes how the coarse grained effective action changes with 
the scale parameter k, relating this scale dependence to the second functional derivative of 
Tfc and the scale dependence of the IR regulator function Rk. The boundary conditions for 
the flow equation are such that the flow relates the microscopic action S = lim^oo Tk with 
the corresponding macroscopic effective action T = lim^o Tfc, the generating functional of 
1PI Green functions. 

To be more explicit, we follow the 'effective average action' approach as advocated in |5j 
and consider the flow equation 

I^H^ 1 * 1 **)" 1 ^}- (2l) 

Here, $ denotes bosonic fields and t — In k the logarithmic scale parameter. The length scale 
k~ l can be interpreted as a coarse-graining scale |J. The right hand side of eq. ( |2.1|) contains 
the regulator function and the second functional derivative of the effective action with 
respect to the fields. The trace denotes a summation over all indices and integration over all 
momenta. The above flow interpolates between the classical and quantum effective action 
due to some properties of the regulator functions Rk (see Sect. |V1B| ). It is important to 
realise that the integrand of the flow equation fl2.1|) , as a function of momenta q, is peaked 
about q 2 « k 2 , and suppressed elsewhere. Consequently, at each infinitesimal integration 
step k — >• k — Ak only a narrow window of momentum modes contribute to the change 
of Tfc — ► Tfc-AAi- In particular, modes with momenta q ^> k do no longer contribute to 
the running at the scale k. It is this property which justifies the interpretation of Tk as a 
coarse-grained effective action with modes q^> k already integrated out. 



For gauge theories, the flow equation (|2.1|) has to be accompanied by a modified Ward iden- 
tity which has to be satisfied at each scale k. Such a requirement is necessary to guarantee 
that the physical Green functions obtained for k —>■ obey the usual Ward identity fl9|-|2~2" 



Here, we use the background field formalism, as employed in || in the context of the effective 
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average action with a covariant gauge fixing (Landau gauge). 

The flow equation couples the infinite number of operators describing an effective action 
with its second functional derivate. In order to solve (|2.1|) , one has to truncate Tk to some 
finite number of operators relevant for the problem under investigation. Some systematic 
expansions for the flow equations are known. Apart from a weak coupling expansion, which 
is known to reproduce the standard perturbative loop expansion, one can use expansions in 
powers of the fields, derivative expansions, or combinations thereof. These latter expansions 
have the advantage of not being necessarily restricted to a small coupling regime. A discus- 



sion on the use of a derivative expansion in Wilsonian RG is presented in p3 |. 

We now turn to our Ansatz for the Abelian Higgs model. The most important information 
regarding the phase structure of the model is encoded in the effective potential (or coarse 
grained free energy) U k , from which all further thermodynamical quantities are derived. 
Equally important is the wave function renormalisation factor of the gauge fields Zp, which 
encodes the non-trivial running of the Abelian charge. In turn, the wave function renormal- 
isation factor for the scalar fields is less important because the scalar field anomalous 
dimension remains small in the type-I region of phase transitions. Hence, we approximate 
the effective action I\. to leading order(s) in a derivative expansion through the following 
operators 

r k [<f>, A] = J d d x (u k {p) + ^ZwF^F^ + Z^D^AjipyD^Ajipl (2.2) 

where p = (p*<p, F^ u = d^A v — d v A^ is the field strength of the electromagnetic field, and 
Dp denotes the covariant derivative — ieA^. 

In principle, the flow equation can be used directly (starting with initial parameters of the 
4c? theory at T = 0) to compute the corresponding critical potential at finite temperature 
within the imaginary time formalism, or, like in f24j, using a real-time formulation of the 



Wilsonian RG ||25|| . Our strategy in the present case is slightly simpler. We are interested 
in the region of parameter space where the 4d couplings are small enough to allow a per- 
turbative integrating-out of the super-heavy and heavy modes, i.e. the non-zero Matsubara 
modes for all the fields and the Debye mode. In this case, we can rely on the dimensional 



reduction scenario and employ the results of | 2~6[] , who computed the initial conditions per 



turbatively. The result is then a purely three-dimensional theory for the remaining light 
degrees of freedom, whose infra-red behaviour is then studied applying the above Wilso- 
nian renormalisation group. In the sequel, we will therefore need the flow equations for Uj. 
and e 2 (k) in 3d. At the scale of dimensional reduction, that is the starting ultra-violet (UV) 
scale A of the 3c?-flow, we normalize the wave function factors to one, and the initial effective 
potential U\ is obtained from dimensional reduction. 



B. Cross-over of the gauge coupling 

We now consider the case d = 3, and discuss the flow for the Abelian coupling. A main 
feature of the Abelian Higgs theory in 3d is that the Abelian charge scales in a non-trivial 
manner with the coarse graining scale k. The dimensionless Abelian charge in 3c? is defined 

as 
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„ 2 ,m _ 4W _ 5§(*0 



and its scale dependence is related to the gauge field anomalous dimension rjp = —d t In Zp(k) 
(here a function of k and the fields) through || 

d ~§ = -4 (1 - Vf) ■ (2-4) 

The first term in Q2.4j) comes from the intrinsic dimension of the charge squared (propor- 
tional to k), while the second term proportional to the gauge field anomalous dimension 
accounts for the non-trivial running of the coupling. The flow (|2.4| ) has always the (trivial) 
Gaussian fixed point given by e 2 = 0. In addition, one might encounter further non-trivial 
fixed points which are given implicitly through the solutions of rjp = 1. 

Both the scalar and the gauge field anomalous dimensions and np are perturbatively 
small near the Gaussian fixed point, i.e. and \r) F \ <C 1. This holds true at the initial 
scale for k = A in the effective 3d running to be specified later. It follows that the running of 
the dimensionful Abelian charge is negligible near the Gaussian fixed point, e 2 (k) « e 2 (A). 
Here, the dimensionless coupling scales as e 2 (k) ~ A/k. In this regime it is expected that 
standard perturbation theory gives a reliable estimate of the effective potential in this region 
of the parameter space |L3] . 



However, for rjp < 1 the Gaussian fixed point is IR unstable, as follows directly from 
Therefore, when approaching the infra-red, the dimensionless Abelian charge will unavoid- 
ably grow large, scaling away from the Gaussian fixed point. In particular, it can enter 
into a region where rjp(e 2 ) is no longer <C 1. When a non-trivial fixed point is approached, 
i.e. rjp fa 1, the suppression factor (1 — rjp) in ( |2.4| ) becomes important. A strong linear run- 
ning of e 2 ~ k (the IR region is effectively 3d) will ultimately set in as soon as the deviation 
from the Gaussian fixed point becomes sizeable P, pH| , p7[ . In this regime, we expect some 
quantitative modifications of the predictions by perturbation theory due to the non-trivial 
running of the Abelian charge. 



C. Abelian fixed point 

The anomalous dimension r]p has been calculated in ||. It is, in general, a complicated 
function of the gauge coupling, the fields, and the further parameter describing the effective 
action in a given approximation, like the coarse grained potential (cf. (113) of |§). However, 
rjp is proportional to e| itself, and we write it as 

Vf(P) = -4-, ■ (2-5) 

Given the anomalous dimension, ( |2.5| ) provides a definition of e 2 (p). Our current under- 
standing of the IR behaviour of the gauge sector hinges on the precise properties of rjp(p), 
and hence of e 2 (p). 

Let us recall a few cases where e 2 is approximately known. First, within standard pertur- 
bation theory, the dimensionful gauge coupling e 2 = e\k = const throughout. Within our 
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formalism, the 'no-running' corresponds to the limit e 2 — > oo. In this limit, the effective 
fixed point is independent of the fields and we can expect to be close to the results from 
perturbation theory, as long as additional effects due to the scalar anomalous dimension can 
be neglected.0 

Second, consider the large -N limit of the £/(l)-Higgs model, where N denotes the number 
of complex scalar fields. In this limit, the flow (|2.4j) is dominated by the contributions of the 
Goldstone modes. They overwhelm those due to the radial mode. Therefore, e 2 becomes 

9 6-7T 2 

<S « -Jjr (2.6) 

close to the minimum of the effective potential.^ In particular, ( |2.6j ) does no longer depend 
on the quartic scalar coupling or the location of the v.e.v. because the massive (radial) mode 
is suppressed. Extrapolating ( j2.6| ) down to the physically relevant case N = 1 corresponds 
to replacing the radial mode by a massless one. This yields e 2 ~ 6tt 2 in accordance with 
the leading order result from the e-expansion. This value serves as a reference value for our 
subsequent considerations. 

Third, we recall the findings of || and ||10|| , where the function rjp has been studied nu- 
merically for different iV within a local polynomial approximation of the flow about the 
non-trivial minimum at p = po (up to ~ <f) s ). It was found that the implicit solutions to 
rjF{el) = 1 for small N (in particular N — 1) can deviate considerably from the large -N 
extrapolation 67r 2 . This deviation is due to the decoupling effects of the massive mode. 
Still, the qualitative form of ( |2.5|) , where the function e 2 is replaced by an effective field- 
independent fixed point, remains a good approximation to (|2.5| ). This simplified picture 
persists if the field derivatives <91n[e 2 (p)]/<9p remain small within the non-convex region of 



the effective potential (see also the discussion in Sect. |VE|) . This implies that the threshold 
effects of the radial mode for N — 1 act on ( |2.4j ) as varying the number of scalar fields in 

m. 



Hence, the qualitative structure of the flow (|2.4|), to leading order, is determined by ( p.5| ) 
with e 2 given by some number, e.g. the appropriate effective fixed point. For the present 
purpose it is sufficient to study the flow Q2.4p with e 2 as a free parameter. The properties 
of the first-order phase transition depend on the size of e 2 . However, as we shall see in 
detail below, the dependence turns out to be very small for large e 2 : this part of the phase 
diagram can be studied without having a complete understanding of the underlying fixed 
point structure. In turn, we find a strong dependence within regions where the effective 
fixed point is small. For this case, a more refined analysis is required in order to provide 
more reliable predictions. 



3 In the region where A>e 2 (e.g. strongly type-II superconductors), the critical behaviour of the 
limit e 2 — > oo corresponds to an effective scalar theory which belongs to a different universality 
class than the 0(2N) scalar theory obtained for e 2 — > ||. 

4 In ([O]), '~' means equality up to a regulator scheme dependent coefficient of 0(1). 
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D. Cross-over scale 



Within the remaining part of the article we approximate the anomalous dimension as de- 
scribed above. Hence, the eqs. (|2.4|) and (|2~5|) are easily solved by 



A 



We note the appearance of a characteristic cross- over scale 

Aej{A) 

" ~ el-el(A) ■ (2 ' 8) 

It describes the cross-over between the Gaussian and the Abelian fixed point, and depends on 
the initial conditions. For k > k CT the running is very slow and dominated by the Gaussian 
fixed point, e\(k) ~ const . This corresponds also to the limit e 2 — > oo. On the other hand, 
for k < k cr the running becomes strongly linear and the Abelian fixed point governs the 
scale dependence, e|(fc) ~ k. The question as to how strong the first-order phase transition 
is affected by this cross-over depends on whether the cross-over scale is much larger (strong 
effect) or much smaller (weak effect) than the typical scales of the transition (see Sect. |V D| ). 



e 2 = e 2 = el e 2 = oo 
I ►> 1 "*~* 1 

k cr > k cr < 

FIG. 1. The relation between the sign of the cross-over scale k cr and the running of 
the gauge coupling (arrows indicate the direction of the e 2 flow as k — > 0). 



The cross-over scale turns negative if the initial value e§(A) is too big. This simply means 
that the flow would never be dominated by the Gaussian fixed point (see Fig. |1]) in the first 
place (no cross-over). Although this case is interesting in its own right, this region will not 
be discussed any further. 



E. The running potential 

We now turn to the flow equation for the effective potential, which can be obtained from 
the flow equation ( |2.1| ) using the Ansatz given by (|2.2p. The resulting flow equation is a 
second-order non-linear partial differential equation. It has been derived originally in 
and reads in 3d 

£ L UM = (w _ 1} a (<M) + tl ( Ml ) + 2 el (?mi) ( , 9) 

for the case of N complex scalar fields. Similar flow equations are obtained for the wave 
function factors and Zp, and thus for the anomalous dimensions 77^ = —dt\nZ 9 and 
7] F = —d t \nZ F . Here, £q(uj) denotes a scheme dependent threshold function defined as 
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<(«) = - («M. + »)/. dy [v{ i r '^l ]n+1 . (2.10) 







These functions have a pole at some uj < and vanish for large arguments. The function 
r(q 2 /k 2 ) is related to the regulator function R k introduced in (Of) through 



R k {q 2 ) = Zq 2 r(q 2 /k 2 ), (2.11) 

where Z denotes either the scalar or gauge field wave function renormalisation.^ 

We can distinguish three different contributions to the running of the potential (|2.9| ) which 
are, from the left to the right, related to the massless scalar, massive scalar, and the gauge 
field fluctuations, respectively. Not all the three of them are of the same order of magnitude, 
though. Indeed, it was already noted that the gauge field fluctuations dominate ( |2.9j ) if 
the quartic scalar coupling A is much smaller than the gauge coupling squared, A/e 2 <C 1. 
This is the case for the physically relevant initial conditions, that is, for the starting point 
of the flow equation ( |2.9|) . Therefore, we can make a further approximation and neglect the 
contributions from the scalar field fluctuations compared to those from the gauge field. The 
flow equation thus takes the form 

Integrating the approximated flow equation allows to control self-consistently whether the 
effects from scalar fluctuations remain negligible or not. It suffices to evaluate the right-hand 
side of ( |2.9|) with U k from the solution of Q2.12 ) to compare the contribution of the neglected 



terms to the running of, say, U k with the leading contributions. It is well known that the 
scalar fluctuations are important for the inner part of the effective potential which becomes 



convex in the limit k — > [ 29| . Therefore it is to be expected that this approximation 



becomes unreliable, within the non-convex part of the potential, at some scale k^ at . 



The solution to fl2.12| ) is the first step of a systematic iteration to compute the solution 



to (|2.9|). The next step would be to replace U k on the r.h.s. of (f2.9|) by the solution to 



fl2.12j ). Proceeding to the next iteration step the scalar fluctuations are eventually taken 



into account. Solving (2^) with U k on the right-hand side replaced by the explicit solution 
of ( [2.12D is much easier than solving (|2.9| ) directly, because the former becomes an ordinary 
differential equation, while the later is a partial one. This procedure can be interpreted as 
an expansion in terms of scalar loops around the gauge field sector. We will mainly use 
the first step in the sequel. In order to estimate the integrated contribution of the scalar 
fluctuations, we will in addition discuss the solution of ( j2.9|) with U k on the right-hand side 
replaced by Ua (see Appendix |C|). 



5 A more detailed discussion of both R k and the dimensionless functions r(q 2 /k 2 ) is postponed 



until Sect. VI B 



S 



F. The coarse-grained free energy 



The coarse grained free energy obtains as the solution to the coupled set of flow equations 
(|2.4|) and ( |2.9| ). In the present case, a solution can be written as 



U k {p) = U A (p)+A k (p) . (2.13a) 

Here, the term A(p) stems from integrating out the 3d fluctuations between the scales A 
and k. With e\[k) from C pTTf) and dllk/dk from fl2.12j ), it reads 

a / -\ 1 [ A r t f, r'{y)y^¥{l + k/k cv ) 

A k{p) = tt^ dk dy - - + const. 2.13b 

2n z Jk Jo y k(l + r) (1 + k/k cr ) + 2ejp 



The constant is fixed by requiring that Afe(O) = 0. In Eq. (|2.13b|) we see that the term 



resulting form integrating-out 3c? effective modes depends on the RS through the regulator 
function r(y) and its first derivative. (Explicit expressions are given in the Appendix |I|). 

The above expressions are enough to study all properties of the phase transitions as functions 
of the parameters of the potential U\. 

We are aiming to use an initial condition at k = A obtained from perturbation theory in 4c!. 
This requires that the parametrisation of the 3d potential U\ is such that the matching 
equates the right parameters. In the universal limit A — ► oo, the effective mass term 
contained via U A is renormalised to U A — > U\ — C\p. For a sharp cut-off, we find explicitly 



C A (e) = & (Ak cr - fc c 2 r ln(A/A ) ) . (2.14) 

For finite A, this corresponds to a finite renormalisation of the parameters of the theory, i.e. 
the mass term, or, equivalently, a finite shift of the v.e.v. at the matching scale.f] This finite 
renormalisation has its origin simply in the way how the flow equation integrates-out the 
3d momentum scales. Only after this transformation it will be appropriate to identify the 
potential U\ at the scale of dimensional reduction with the renormalised effective potential 
obtained from a perturbative calculation. 

III. THERMAL INITIAL CONDITIONS 

We now specify in concrete terms the initial conditions for the effective 3d theory. The task 
is to relate the 3c? renormalised parameters of the effective potential to those of the T = 
4c? theory. The initial conditions for the 3d running are the potential U\(p) and the gauge 
coupling (A). The effective perturbative 3c? Lagrangean has been derived in |26[. We start 
with the 4c! effective action, 

A]=JdS\\ F, V F, U + (V,<t>y(V,cf>) - ^ 0V + ^ (0V) 2 } , (3.1) 



6 This shift corresponds to the finite renormalisation as employed in [12| 



9 



where is a single component 4c? complex scalar field. The mass parameter mn entering 
( |3.1|) denotes the T = Higgs boson mass. It is related to the other zero temperature 
parameters of the theory by 

A rriu 

# = 4 < 3 - 2 ' 

with Mw the photon mass. In the phase with spontaneous symmetry breaking, > 0, we 
have (<f)*<f)) = v 2 /2 = M^/2e 2 . The effective action for the 3d theory obtains as 

T A [<p, A]= J d 3 x {^FijFij + (Vitfpitp) + V A (p)} , (3.3a) 

y A (p)=m^V + ^ ( ^ )2 (3 . 3b) 

where ip is the static component of and i,j the spatial components of /i, v. The electric 
component of the gauge field has been fully integrated out because it acquires a thermal 
(Debye) mass m^. The effects of the fluctuation of this mode are suppressed by inverse 
powers of T as oc T, like the nonstatic modes. Following | 26j , the matching conditions 
read to 1-loop accuracy 

e 2 3 (A) = e 2 T (3.4a) 

(4 \ 4 rp2 
A + ^—z ] T - =-— (3.4b) 
4vr 2 j 4vr m D (A) v ' 

mt(A) = Qe 2 + iA) T 2 - l -ml - ^Tm D (A) (3.4c) 
m 2 ) (A) = ie 2 T 2 . (3.4d) 

Using the above, and taking into account the finite renormalisation ( [2.14| ) as explained in 
Sect. [II F| , the renormalised effective initial potential U\(p) entering ( j2.13a| ) can be expressed 
in terms of the T = parameters and Q3.4|) as 



with 



U\(p) = -m 2 l p+ ~ A R p 2 (3.5a) 



-4(A) = \ < - ( j + - 6 - j^J + cm , (3.5b) 

MA)=(A+^-^)r, (3.5c) 

and the dimensionless renormalised quartic coupling reads Ar = Ar/A. The renormalised 
v.e.v. at the scale of dimensional reduction follows as 

Pr(A) = m R (A)/A R (A) . (3.5d) 
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All the 3d parameters are now denned at the reduction scale A, which is on dimensional 
grounds linearly related to the temperature, 



A = £T. 



(3.6) 



Using eqs. (|2.S| ), fl3.4a|) and (|3T6| ) it follows, that the cross-over scale k CT is also related to T 

as 

£e 2 



iel-e 2 



T . 



(3.7) 



Let us finally comment on the matching parameter £. On one hand, £ has to be smaller 
than 2ir, because elsewise the assumption that all heavy modes have been integrated out 
can no longer be maintained. On the other hand, a too small value for £, say £ < 1, would 
tend to neglect contributions from modes roughly within the window ~ 2nT and ~ T. For 
the problem under consideration £ ~ 1 turns out to be a good choice. This choice shall be 
adopted throughout. Our results do depend very little on a variation of this matching scale 
(see also the comment in Sect. [V D| below) . 



IV. THE PHASE DIAGRAM AT FINITE TEMPERATURE 

We have now all the ingredients to study in detail the phase diagram and the phase transition 
of scalar electrodynamics. In this section, we discuss the main characteristics of the phase 
diagram as well as some properties of the critical line. The following section collects our 
results for the thermodynamical quantities related to the first-order phase transition and a 
discussion of the characteristic scales of the problem. 



T 

M 

4 



3 



2 



1 











SYM 






SSB 



20 40 60 80 100 120 

mn[GeV] 



FIG. 2. The phase diagram in the (T, mn)-plane. 
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A. The phase diagram 



The 'phases' of scalar electrodynamics are distinguished by the location of the global mini- 
mum of the effective potential. Above the critical temperature, the ground state corresponds 
to vanishing field po = 0, that is, to the symmetric phase (SYM). Below the critical tempera- 
ture, the ground state corresponds to po 7^ 0, the phase with spontaneous symmetry breaking 
(SSB).[| The corresponding phase diagram in the (T, mn)-plane is displayed in Fig. ||]. The 
phase transition between these two phases is first order for small \$/e\, that is for small 
values of the Higgs field mass. In the context of superconductivity this region corresponds 
to the strongly type-I systems. For very large Higgs field mass, the phase transition turns 
second or higher order .[] 




0.002 0.004 0.006 0.008 0.01 " 7 0.002 0.004 0.006 0.008 0.01 

p/T 

FIG. 3. The coarse grained free energy as a function of the scale parameter and the temperature 
(mn = 60 GeV). The full line corresponds to T = T c and k Stable- Left panel: T = T c , for 
different scales k. Right panel: k — ► Actable > for different temperatures around T c . 

In Fig. [3], we have displayed the coarse grained free energy within the type-I region of param- 
eters for rrtH = 60 GeV for different scales and temperatures. At the critical temperature 
(left panel), it is realised that a barrier is building up for decreasing scale k, eventually 
creating a second minima at vanishing field. The minima are degenerate in the infra-red 
limit k « Actable (which corresponds roughly to k — > in the present approximation). Notice 



7 It is sensible to speak of two distinct phases only for N > 1 complex scalar fields. For N = 1, the 
symmetry is never broken in the strict sense. However, we will stick to the usual - albeit slightly 
incorrect - terminology even for N = 1. 

8 The strongly type-II region has been studied using flow equations within a local polynomial 
approximation in Q. See also [12]. 
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that the flattening of the inner part of the potential is not observed because the scalar fluc- 
tuations have been neglected at the present state. Rather, the effective potential reaches the 
degenerate shape already at some scale Arabic, which should be larger than the scale where 
the flattening sets m.f\ 

The temperature dependence of the coarse grained free energy at k A; sta bie is shown in the 
right panel. The metastability range AT = T s — Tj, between the barrier temperature TJ,, 
where the potential develops a second minimum at the origin (lowest dashed curve) and the 
spinodal temperature T s , where the asymmetric minimum disappears (upper- most dashed 
curve), is very small. 

The physical quantities that characterise a first-order phase transition (except the metasta- 
bility range) are defined at the critical temperature T c , when the potential has two degenerate 
minima, the trivial one at p = and a non-trivial one at p = po ^ 0. The critical line of the 
phase diagram as depicted in Fig. El is obtained solving the criticality conditions 







dU k 



dp 



(4.1a) 



P=P0 



£4(0) = U k (po) 



(4.1b) 



Here we kept k arbitrary though strictly only for k = are these conditions required physi- 
cally. They establish a relationship between the parameters of the theory, and thereby define 
the critical line between the symmetric and the SSB phase in Fig. It is helpful to rewrite 
the conditions (14.11) into 



F 1 (p/T) = A R 



mi 



F, {P/T) = 2^. 
The functions iq and F 2 are related to the fluctuation integral through 



2 rr 



A(x) - xA'(x) 



with 



r 

F 2 (x) = — \2A(x)-xA'(x) 
x 1 



A(p/T) = A(p)/T 3 



(4.2a) 
(4.2b) 

(4.3a) 
(4.3b) 

(4.3c) 



The first condition determines the ratio x = p/T c of the discontinuity to critical temperature 
in dependence on the 4c? parameters as given through Aft(e, A) from ( |3.5c| ). The second one 



relates the solution of (|4.2a ) to the ratio of the Higgs boson mass to critical temperature 
and (|3.5b|) , and eventually to the critical temperature and the discontinuity itself. 

Explicit expressions for the scale-dependent effective potential and the function A(p) are 
given in the Appendix M. 



9 'A quantitative discussion of these scales is given in Sect. VP below. 
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B. Endpoint of the critical line 



Some simple properties of the solutions to (|4.2| ) can be deduced directly from the functions 
Fx t 2- For x > 0, these functions [with from ( [B.5| )1 are positive, finite, monotonically 
decreasing and vanishing for x — > oo. They reach their respective maxima at x — 0, with 
(for k = 0) 

F 1 (0) = \e\e 2 , (4.4a) 
= (l - T^~, In (4) ) ■ (4.4b) 



71-2 £ e * — e2 V £ e * — e 



The renormalised 3<i quartic coupling Ar, as given by (|3.5cj ) and fixed through the 4d pa- 



rameters of the theory, is positive in the domain under consideration. Given the monotony 
property of Fi it follows that a solution to (|4.2a|) is unique (if it exists). There exists 



no solution for too large values of Ar. Its largest possible value corresponds to vanishing 
v.e.v. i.e. to x = 0. Using eqs. (|3.2|) and (|3.5c| ) gives an upper bound on the scalar mass for 



the phase transition being first order. It reads 

4 - ^ • (45) 

For any finite value of e\ ( |4.5|) predicts an upper limit for the mass of the Higgs particle. 
This is an immediate consequence of the existence of an effective fixed point for the running 
gauge coupling ( |2.4| ). Indeed, as the limit e\ — > oo corresponds to perturbation theory we 
recover the standard perturbative prediction of a first-order phase transition for all Higgs 
boson mass. This endpoint is usually interpreted as the tri-critical point of the model, 
above which the phase transition turns from a first-order transition to a second-order one. 
However, the endpoint of the first-order transition line is within the domain of validity of 
the present computation only for sufficiently small values of e^.Q For larger values of the 
Abelian fixed point, we expect that the precise location of the endpoint is also determined 
by the scalar field fluctuations. 

In the opposite case, the smallest possible value for Ar corresponds to x — > oo, thus to 
Ar = 0. This gives a lower bound on the mass of the Higgs particle according to 



2 



^£<^h (46) 

4tt < M*, ■ l4 - DJ 

For My/ = 80.6 GeV and e = 0.3 the bound is satisfied at about ttt-h ~ 16 GeV. This bound 
stems entirely from the initial conditions employed. This indicates that the dimensional 
reduction scenario is no longer appropriate for small mn- In the present work, we are also 
not interested in the region of parameter space where the Coleman- Weinberg mechanism 
already takes place within the original 4d theory, which happens at even smaller values for 
777 H (typically for A/e 4 at about 3/87T 2 or smaller p7[). 



10 The endpoint presented in Fig. ^ corresponds to el ~ 5 
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V. THERMODYNAMICS OF THE FIRST-ORDER PHASE TRANSITION 

Here we present our results for the coarse grained free energy and related physical quantities 
close to the critical temperature of the first-order phase transition as a function of the 
effective Abelian fixed point. The initial conditions are specified through the gauge coupling 
at vanishing temperature e = 0.3 and the photon mass Mw = 80.6 GeV. The ratio A/e 2 of 
the Ad couplings ranges between 0.06 — 0.75 for a Higgs field mass between 20 — 70 GeV. 
Relevant information is given by the critical temperature T c , the discontinuity at the phase 
transition p , the latent heat L and the surface tension cr^ We compare our findings to 
perturbation theory, and to lattice simulations (for the critical temperature). All our results 
are obtained as functions of the effective fixed point of the Abelian charge. Due to the 
approximations performed, they depend also on the regularisation scheme. We use a sharp 
cut-off regulator throughout the present section. The regularisation scheme dependence is 
discussed in the following section. 

A. Discontinuity and critical temperature 

We begin with the discontinuity and the critical temperature, which follow directly from 
solving the criticality conditions ( f4.2p . 



m H = 30 GeV 



T c /ra H 



m H = 50 GeV 



m H = 70 GeV 



2 5 10 20 50 100 

FIG. 4. The critical temperature as a function of the Abelian fixed point. 
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A comment concerning the dimensions is in order: U, a, L and p will be given in 3d units, unless 
otherwise stated. Their 4d counterparts are simply obtained by multiplying with T. 
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The critical temperature as a function of the Abelian fixed point is given in Fig. f| for 
m H = 30, 50, 70 GeV. It turns out that T c is rather insensitive against e*. We observe an 
effect of a few percent only for very small values of el (see also Fig. [15]). This is not a feature 



of the Higgs mass being relatively small, as similar results are obtained for all m^. 

Before continuing, let us briefly compare our findings for the critical temperature to existing 
lattice data. Lattice results have been reported for e = |, mw = 80.6 GeV and mn = 30 
GeV for the non-compact [/(l)-Higgs model in [I7|, and for the compact one in [18]. The 



result reported in |L7j is T c = 131.18 GeV for a finite lattice spacing. The continuum limit 



gives the slightly lower value T c = 130.86 GeV [18]. This is consistent with T c = 131.28 



GeV, the result for the compact case |I8| . Here, for el = 6rr 2 , we find T c = 128.11 GeV. 



As follows from Fig. [|, the critical temperature is essentially independent of the effective 
Abelian fixed point. The perturbative value is T c = 132.64 GeV [17||. These results are in 
good numerical agreement. 

We now turn to the discussion of the discontinuity. In Fig. |5| we compare the logarithm of the 
v.e.v. in 4d units (normalised to the v.e.v. at T = 0) at different scales. The renormalisation 
of po between the T = and the k = A lines results from the integration of the heavy and 
super- heavy modes, given by fl3.5d| ). The scale k vcv is defined as the scale where the running 
of the potential minimum stops. This scale is related to the scale /cm, where the photon mass 
in the SSB regime is becoming larger than the coarse graining scale, and thus decouples. 
Indeed, in the present approximation, the flow equation for the potential minimum reads 



dp 1 e 2 {k) 



il(M 2 (k)/k 2 ) . (5.1) 



dk it 2 X(k) 

Here, X(k) = U')!(p (k)) denotes the quartic coupling at the minimum, and 

M 2 {k) = 2e 2 {k)p {k) (5.2) 

the photon mass squared. The running of the v.e.v. decouples at k « k vcv , which happens as 
soon as the 3d photon mass M is sufficiently larger than the scale k (roughly at M 2 /k 2 « 10) 
such that the threshold function in (|5.1| ) suppresses any further renormalisation (k vev /T is 
displayed in Fig. |T^). 

From Fig. [| we conclude that the main part of the actual running of the potential minimum 
comes from integrating-out the 3d fluctuations, as can be inferred from the wide separation 
of the k = A and the k ~ k vev lines as opposed to the comparatively narrow separation of 
the T = and the k = A lines. 

Fig. ^ shows the v.e.v. po as a function of the Higgs field mass and the Abelian fixed point. 
The shaded region covers the region 2 < e* < oo for the Abelian fixed point value. For 
small mH, the effect is clearly negligible. With increasing mn, however, the influence of the 
running gauge coupling is increasing drastically, leading to a strong weakening of the phase 



transition (see also Fig. |15|) 
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Ad: T = 



3d: k = A 




20 30 40 50 60 70 
mn[GeV] 

FIG. 5. The size of the v.e.v. po(k) at T = 0, k = £T and at k ~ &Vev (in 
units of the 4d v.e.v. at T = 0). 
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FIG. 6. The 3d v.e.v. po as a function of the Abelian fixed point. 

Finally we compare in Fig. [7] the ratio of the Ad v.e.v. 0o to the initial T = v.e.v. v/\/2 for 
different Abelian fixed point values with the findings of perturbation theory.^ Again, the 
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We thank A. Hebecker for providing his data from IS] for comparison in Figs. [7], ^, [H] and 11 
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shaded region covers the region 2 < e* < oo for the Abelian fixed point. We observe that 
the v.e.v. shows a small dependence on the Abelian fixed point for sufficiently small Higgs 
field mass. For larger values of mn, the v.e.v. approaches the perturbative two-loop result. 
It follows that the v.e.v. is rather stable against effects from the running Abelian charge, 
say a least for el > 20. Only for e 2 ~ 4 the running becomes strong enough to result in a 
significant decrease of po- 




B. The critical potential 



The critical potential is shown in Figs. § and |J] for mn = 38 GeV. Fig. ||| gives the critical 
potential in units of the critical temperature for different values of el as functions of 
We note that for large el > 6n 2 , the shape of the potential is rather insensitive against a 
change in e 2 .. Here, the additional scale dependence induced through the gauge coupling 
is quite small (a few percent). For small values of el, the height of the barrier is reduced 
significantly, up to a factor of 3 at el = 4. The strong scaling of e 2 thus weakens the phase 
transition considerably for small el -C 6tt 2 . Again, the quantitative change depends strongly 
on the value for the effective Abelian fixed point, if el <C Qtt 2 . The non-trivial running of 
e 2 (k) has a stronger effect on the small p region of the potential. Here, the decoupling of 



13 Notice that comparing critical potentials (or other relevant quantities) in units of T for different 
values of e 2 is sensible due to the very weak dependence of T c on the effective fixed point (see 
Figs. | and |l5|) . 
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the gauge field sets in only at smaller scales, which in turn results in a stronger quantitative 
effect due to the running gauge coupling. 




0.02 0.04 0.06 0.08 



PIT 

FIG. 8. The critical potential for different values of the effective fixed point. 




0.05 0.1 0.15 0.2 0.25 



y/2<j>/v 

FIG. 9. The critical potential, comparison with perturbation theory (see text). 
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Fig. gives the critical potential in units of the 4c? v.e.v. v /y/2, and compares the solution 
of (|2.12 ) with those obtained within perturbation theory (PT). Line (a) corresponds to PT 
to order (e 3 , A 3 / 2 ) fl4]| , line (b) to our result with e 2 = 6tt 2 , line (c) to PT at order (e 4 , A 2 ) 
15] and line (d) to our result with e 2 = 4. For e 2 = 6n 2 , the critical potential is situated 



half way between the one- and two-loop perturbative results. For decreasing e 2 , the critical 
potential approaches quickly the two- loop result, and becomes even smaller at about e 2 ~ 4. 
It is interesting to note that a value for e 2 can be found for which the two-loop perturbative 
result is matched perfectly. 



C. Surface tension and latent heat 

The interface tension for a planar interface separating the two degenerate vacua follows from 
as 

a = 2 J o d<p^J Z^Ucritip) • ( 5 -3) 

It is sensible to the actual shape of the critical potential and yields additional information 
regarding the strength of the phase transition. In Fig. [H] the surface tension is shown as a 
function of mn and in comparison with perturbation theory. The shaded region covers the 
results for 2 < < \^Qtc. We again note that the effect of the running coupling is negligible 
for small Higgs boson mass. In contrast to the v.e.v., the surface tension depends rather 
strongly on e 2 already for moderate values of mn- An even stronger running of e 2 would 
lead to a dramatic decrease of the surface tension, up to several orders of magnitude. 




20 30 40 50 60 70 
mn[GeV] 



FIG. 10. The surface tension for various values of the Abelian fixed point 
in comparison with perturbation theory to order (e 3 , A 3//2 ) and (e 4 , A 2 ). 



20 



Finally, we consider the latent heat L, defined at the critical temperature as 



Using eqs. 



L = T f dUjpo) _ dU(0) \ 
\ dT dT ) 

3[) and (EToT) we obtain 



(5.4) 



T=T C 



L = (m| - 2mjQ p + t^rTPq + 3A (Po) - PoA'(p 



(5.5) 



The latent heat is related to the discontinuity and the mass of the scalar particle. Using 
04.21), it can be shown that 



Po Wh 



(5.6) 



which is also known as the Clausius-Clapeyron equation. This relation was shown to be 
fulfilled within an explicit gauge- invariant perturbative calculation |16|]. However, it holds 
not true within standard perturbation theory: the perturbative values for the latent heat 
as found in ]15| are all below the value given through the Clausius-Clapeyron relation ( f5.6|) . 
The deviation varies between a few percent up to 15-20% for mn between 20 GeV and 
70 GeV, and is larger at order (e 4 , A 2 ) than at order (e 3 , A 3 / 2 ). 
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FIG. 11. The latent heat for various values of the Abelian fixed point in 
comparison with perturbation theory to order (e 3 , A 3//2 ) and (e 4 , A 2 ). 
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The latent heat in units of the critical temperature is displayed in Fig. [11] for various values 
of the effective Abelian fixed point, and in comparison with perturbation theory to order 
(e 3 , A 3 / 2 ) and (e 4 , A 2 ). The shaded region covers the interval 2 < e* < \^6n. We again observe 
a sharp decrease for small e* and large Higgs boson mass as in Fig. [?]. It is interesting to 
note that the curve for e* = 4 roughly agrees with the two-loop perturbative result for all 
inn above 30 GeV. This is not the case for the surface tension. Comparing Fig. [TT] with 
Fig. ITTjl we notice that the effect of the running gauge coupling is more pronounced for the 



surface tension, because the entire region for p < po enters ( |5.3D , while the latent heat is 
only affected by po. 



D. Characteristic scales 

We discuss the results obtained so far in terms of the characteristic scales relevant for the 
phase transition. Most of the qualitative (and even quantitative) features can be understood 
once these scales are known. 



In Fig. 12, we have depicted the relevant momentum scales as a function of the Higgs mass. 
The top line at k = A corresponds to the scale of dimensional reduction, that is, the starting 
point of the flow in 3d. The scales k s ,k vev and Stable (full lines) describe characteristics 
of the potential, the scale k cr (dashed lines, for two values of the Abelian fixed point) the 
characteristics of the gauge sector, and /cfl at (dashed-dotted line) the scale where scalar fluc- 
tuations can no longer be neglected within the non-convex part of the potential. All these 
scales are now discussed in detail. 

At k = k s , the origin of the effective potential stabilizes, U'(p = 0) = 0, as the mass term 
squared at vanishing field changes sign. The free energy has two local minima for scales 
below k s . This scale is therefore a good estimate for the scale of discontinuity. In ||, an 
estimate for this scale has been given, based on a local polynomial approximation for the 
potential. Within our conventions, it reads k^ s « 0.18c 4 (T)/Ar(T) for a sharp cut-off, and 
roughly coincides with k s as presented here (k dis /k s ranges between l to 3). 

The scale k ~ k vev indicates when the v.e.v. p is within 1% of its final value, eventually 
reached for k — > 0. However, this is not yet the scale where the critical potential has reached 
a stable shape, which actually happens only at about k « Stable • This results from the fact 
that the effective photon mass squared 2e 2 (k)p (within the non-convex part of the potential) 
is smaller than the photon mass at the minimum in the SSB regime Q5.2Q , and the decoupling 
takes place only at smaller scales. Here, we have obtained Stable comparing the depth of the 
potential U(0) — U(p ) at p with the height of the barrier ?7(p max ) — U(Q), demanding this 
ratio to be below « 5%. At k — k stllUe , the v.e.v. is as close as 0.1% to its final value. 

The cross-over scale k cr characterises the cross-over from the Gaussian to the Abelian fixed 
point. For e 2 = 6n 2 , we see that k CI is about 1-2 orders of magnitude smaller than the scale 
k s , which explains why the running gauge coupling has, in this case, only a small numerical 



Remember that the critical potential at fc s tabie> within the present approximations, is about the 
same as at k = 0, as no substantial running takes place below fc s tabie- 
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effect on the properties of the phase transition. From the fact that the scales k vev and fc s t a bie 
are separated by an order of magnitude {k vev / k stah \ e ~ 5), we can conclude that the run- 
ning of the gauge coupling has a stronger effect on physical observables based on the entire 
effective potential (like the surface tension), than those related only to the v.e.v. (like the 
latent heat). This is quantitatively confirmed by the findings displayed in the Figs. ^ p|, |ICi| 
and [II]. For el = 4, we realize that the corresponding cross-over scale is of the same order 
of magnitude as the scales k s , k vev and Stable -Q This is the region where the running of 
the gauge coupling has a strong quantitative effect on the properties of the phase transition, 
leading to a significant decrease of the strength of the transition. 
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FIG. 12. Characteristic scales (see text). 



Finally, we have also indicated the scale k^ (dashed-dotted line), which is an estimate for 
the scale where the flattening of the inner part of the effective potential sets in. We obtained 
/cfiat from solving k 2 + U' k {p) ~ numerically for k in the non-convex part of the potential, 



15 In Fig. 12, the scales k s , /c vev , Actable and k^t have been obtained for e 2 = 6ir 2 . The correspond- 
ing results for e 2 = 4 do deviate (for larger Higgs mass) only slightly from the curves as presented 
here. This minor difference is of no relevance for the present discussion. 



23 



with Uk the leading order solution for the free energy.^ 

In [I3| an estimate for the ratio of fcfl a t/&stabie has been obtained, based on an investigation 
of the surface tension of the 3d Abelian Higgs model in the universal limit A — > oo. There, 
it was found that fcl at / Stable ~ e 2 /M, with M being the 3d photon mass. The boundary 
^LtAstabie ~ 1 yields the relation fcfl at ~ (e 2 T/2po) 1//4 fc s tabie, which, using the data for 
^stable as in Fig. [12], coincides within a few percent with the line for fc flat as obtained above. 
Corrections to the universal limit can be expanded as a series in M 2 /A 2 ||13|| . In the present 
case, we start at a finite scale A = £T, but the smallness of M 2 /A 2 (ranging from 0.2 to 
0.001 for 20 GeV < mn < 70GeV) is responsible for the small corrections with respect to 
the universal limit A — ► oo. Being close to the universal limit of the effective 3d theory also 
explains why the dependence on the matching parameter £ is rather small. 



We now come back to the discussion of ka at from Fig. [12], which, by definition, sets the scale 
below which the scalar fluctuations trigger the flattening within the non-convex part of the 
potential, and hence the scale below which these fluctuations should no longer be neglected. 
First notice, that the scale of discontinuity k s is bigger than k^ at by an order of magnitude. 
We can thus expect that the scale of discontinuity is only weakly affected by the scalar 
fluctuations. Also, k vev > ka at by a factor of ~ 5. Finally, for small Higgs field mass, fcfl at is 
also smaller than fc s t a bie- In this region, only small quantitative changes are expected if the 
scalar fluctuations are taken into account. This is no longer the case for large Higgs field 
mass, where k^ at > Stable- However, as these effects concerns mainly the non-convex part of 
the potential, and thus quantities like the surface tension, we can still expect that the latent 
heat and the v.e.v. are only moderately affected. 

These last observations are also relevant for the applicability of Langer's theory of bubble 
nucleation. The concept of an interface tension, as defined in ( |5.3| ), is based on the implicit 
assumption that the scale fc s tabie can indeed be identified. A criterion for this being the case 
is the smallness of the perturbative expansion parameter. From our consideration we can 
conclude that this will become more and more difficult for increasing e 2 T c /2p > 1, that is, 
for very weakly first-order phase transitions.^ 



E. Higher order corrections 



Finally, we comment on the higher order corrections which are expected from operators 
neglected within the present approximation. Clearly, the results presented here are affected 
by the approximations performed, most notably through (i) the derivative expansion, (ii) 
neglecting the scalar field fluctuations as opposed to the gauge field ones, (Hi) approximating 
the infra-red regime of the Abelian charge by an effective fixed point, and (iv) computing 
the initial conditions perturbatively. We discuss these approximations now one by one. 



16 A similar though slightly shifted curve for /cfl at is obtained from solving k 2 + U' k (p) + 2pU'^.(p) ~ 0. 
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The treatment of very weakly first-order transitions based on coarse grained potentials has been 
considered in [31]. 
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(i) The leading order terms of the derivative expansion are known to correctly describe criti- 
cal equations of state and scaling solutions for a variety of 0(A r )-symmetric scalar models in 
3d. Although little is known about the convergence of such an expansion, it appeared that 
the smallness of the anomalous dimensions controls the influence of higher order derivative 
operators in the effective action. Therefore, an a posteriori consistency check for the relia- 
bility of the derivative expansion consists in computing the corresponding scalar and gauge 
field anomalous dimension rj^ and rjp- In the present case, this involves more complicated 
higher order threshold functions (for their definitions and further details, see 0). At the 
scale k ~ Stable; we can compute the scalar anomalous dimension self-consistently from 
the explicit solution for the effective potential, obtained while neglecting 77^. We find that 
1 77^,1 < 0.005 in the interval considered, which is consistent with our initial approximation 
r] v = and justifies the derivative expansion within the scalar sector. For N = 1, the gauge 
field anomalous dimension r]p can be estimated in a similar way. It becomes of order one 
only when the non-trivial fixed point is approached. We find that rjp ranges from 0.03 to 0.4 
within the range of Higgs field masses considered here and for el ~ 6tt 2 . A main difference 
between the scalar and the gauge field sector is that the gauge field anomalous dimension 
grows large (rjp = 1) at a scaling solution. Therefore, one expects that higher order correc- 
tions within a derivative expansion (or the momentum dependence of the gauge coupling) 
can become important at a scaling solution and should not be neglected. In the present case, 
however, the scales relevant for the first-order phase transition have been reached before the 
Abelian charge finally runs into its non-trivial fixed point, that is before T]f = 1. Therefore 
we can expect that the derivative expansion behaves reasonably well even for the gauge field 
sector. 

(ii) In the same way, we can check the validity of neglecting scalar fluctuations within 
the non-convex part of the effective potential. It is found that the self-consistent inclusion 
of scalar fluctuations to leading order results in corrections of the order of a few percent, 
increasing with increasing Higgs field mass (see Appendix £J). This agrees also with the 
discussion of the preceding section, where it was argued that scalar fluctuations should no 
longer be neglected as soon as fcfl at is of the order of Stable- Clearly, the weaker the first-order 
phase transition the more scalar fluctuations will become relevant at the phase transition. 
For a quantitatively more reliable computation of thermodynamical quantities in the weakly 
type-I region, one has to go beyond the present approximation and include scalar fluctu- 
ations. All the present approximations can be improved in a systematic way, as has been 
emphasized earlier. This can be done either along the lines outlined in Sect. |11 E| , or by a 
straightforward numerical integration of the flow equation as in |12j| . 

(Hi) The main uncertainty in the present understanding of the [/(l)-Higgs theory is linked 
to the gauge sector of the theory i.e. the precise infra-red behaviour of the Abelian gauge 
coupling. Here we have effectively parametrised this uncertainty in terms of an Abelian 
fixed point motivated by previous work based on large -N extrapolations and Wilsonian RG 
techniques. A precise determination of the correct fixed point requires the study of the mo- 
mentum and of the field dependence of the Abelian charge. Our approximation assumes that 
the field gradients of the function el(k,p) remain sufficiently small within the non-convex 
part of the potential at scales above k ~ fc s top- In the large -N limit, where this fixed point is 
well understood, the results in the present approximation are in very good agreement with 
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the result found within a fixed dimension computation. 

(iv) The points (i) — (in) concerned the approximations on the level of the flow equation. 
These are the most important ones, because they act back on Tk upon integration of the 
flow. An additional approximation concerns the initial conditions to the flow. Here, they 
have been obtained from the dimensional reduction scenario within a perturbative loop com- 
putation. For the present purposes, it was sufficient to use a 1-loop perturbative matching 
as given in Sect. [TTT], although the 2-loop matching has been reported as well JHJ. These 
higher order effects can be taken into account in principle; in practice, this shall not be 
necessary because their quantitative influence is smaller than the contributions from the 
scalar fluctuations for larger Higgs field mass, which have already been neglected. In any 
case, a small change of the initial condition cannot change the main effect reported here. 
Except for small Higgs field masses, the dominant contributions come from integrating-out 
modes in 3d This follows directly from Fig. |5|, which shows that the main running of the 
v.e.v. takes place below the scale of dimensional reduction. 

Finally, we remark that the quality of a given approximation can also be assessed by study- 
ing the dependence on the coarse graining scheme. This discussion will be the subject of 
the following section. 

VI. SCHEME DEPENDENCE 

All quantitative results present up to now have been obtained for a sharp cut-off regulator. 
It is a straightforward consequence of the Wilsonian renormalisation group approach that 
physical observables obtained from a solution to a Wilsonian flow equation will not depend 
on the precise form of the coarse-graining. Unfortunately, this conclusion holds only if the 
full effective action is computed. On a technical level, this is barely possible, and truncations 
of the effective action have to be employed. It is precisely this truncation that can introduce 
a spurious coarse graining scheme dependence for physical observables. In this section we 
address the question as to what extend the physical observables as obtained in the preceding 
section do (or do not) depend on the precise form of the coarse graining. In doing so, we 
are able to present quantitative 'error bars' related to the scheme dependence. We also 
present evidence for an intimate quantitative link between the scheme dependence and the 
truncations employed. 

A. Scheme dependence vs. truncations 

Consider the case of computing some physical observable from the solution to a (truncated) 
Wilsonian flow. It goes without saying that a strong dependence of this observable on the 
coarse graining employed is not acceptable as it would cast serious doubts on the truncations 
performed so far. With 'strong' we mean 'inducing large quantitative', or even 'qualitative' 
changes. On the other hand, a weak scheme dependence of certain physical observables is 
a sign for the viability of the approximation employed. In fact, if we were able to solve the 
flow equations without truncating the effective action T^, the final result in the physical 
limit k — > 0, which is by construction nothing else but the full quantum effective action T, 
should not depend on the details of the particular coarse graining employed. There is little 
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hope for this holding true for any truncation of the effective action r& as any truncation 
necessarily neglects infinitely many operators. 

The coarse graining procedure is implemented through the momentum-dependent operator 
Rk(q 2 ). It couples to all the operators present in Tk in a well-defined way, that is, according to 
the flow equation ( j2.ll ). Replacing a coarse-graining by another coarse graining implies that 
the effective coupling of Rk(q 2 ) to the operators contained in the effective action changes 
accordingly. A truncation of the effective action amounts to neglecting infinitely many 
operators to which the coarse graining, in principle, is sensitive. Therefore, studying the 
scheme dependence will probe whether some relevant operators (for the problem under 
investigation) have been neglected, or not. In this light, the indirect feed-back of some 
relevant operators should manifest itself through some strong eigenmode with respect to a 
change of the coarse graining procedure. 

Although these arguments, as presented so far, are of a purely qualitative nature, we will 
show in the sequel that they can indeed be given a quantitative meaning. 



B. Coarse grainings 

Before studying in detail the scheme dependence of our results, we will briefly review the 
main requirements for a viable coarse graining procedure. There are basically three key 
points to be considered. The first one concerns the possible zero-modes of the propagators, 
which typically cause strong infra-red problems within perturbative loop expansions in d < 4 
dimensions. These are properly regularised, if 

lim R k (q 2 ) > (6.1a) 

holds true. This way, the effective inverse propagator for a massless mode reads q 2 + Rk{q 2 ), 
and has a well-defined infra-red limit. The second point concerns the infra-red limit of the 
effective action 1^, which should coincide with the usual effective action for k — > 0. This is 
the case, if 

lim R k {q 2 ) = . (6.1b) 

k^O 

Finally, we have to make sure that the correct initial effective action in the ultra-violet limit 
is approached which is guaranteed by 

lim Rk(q 2 ) — > oo . (6.1c) 

k^oo 



Any function Rk{q 2 ) with the above properties can be considered as a coarse graining |5,13 
It is convenient to re-write Rf. in terms of dimensionless functions r(q 2 /k 2 ) as 

R k (q) = Zq 2 r(q 2 /k 2 ) , (6.2) 

where Z corresponds to a possible wave-function renormalisation (Z^ = 1 in our approxi- 
mation) . 
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Let us introduce two classes of regulator functions which are commonly used in the litera- 
ture. The first one is a class of power-like regularisation schemes given by the coarse-graining 
function 

r P (y) = y~ n , (6.3) 

and y = q 2 /k 2 . The particular case n — 1 corresponds to a mass-like regulator Rk ~ k 2 , 
and n = 2 to a quartic regulator Rk ~ k 4 /q 2 . These algebraic regulators are often used 
because the related threshold functions can be computed analytically. On the other hand, 
these regulators decay only algebraically for large momenta, which can in principle lead to 
an insufficiency in the integrating-out of the hard UV modes. 

A second convenient class of regulators consists of exponential ones, parametrised as 

r M = ( l n , 1 , (6-4) 

exp(cy n ) — 1 

where c is a constant. The exponential regulator with n — c — 1 has been used previously 
in various numerical investigations P, |I0| , |2"T|| . The suppression of large momentum modes 
q 2 ^> k 2 to the flow is now exponential and thus much stronger than in the case of algebraic 
regulators. It is expected that this property is at the basis for a good convergence of 
approximate solutions. 

Both classes of regulator functions depend on the parameter n, with 1 < n < oo. In the 
limit n — > oo, they both approach to what is known as the sharp cut-off regulator, given by 

r » fe) = w^t> - 1 ■ (6 - 5) 

We will now consider the dependence of certain physical observables on particular choices 
of these regulators. 



C. Tri-critical point and large -N limit 



We have given an estimate for the endpoint of the critical line in ( |4.5| ). Its mere existence 
is closely linked to the presence of an Abelian fixed point, although it will be within the 
domain of validity only for small values of the latter. Both functions Fi and F 2 depend 
explicitly on the RS, and so does the solution to eqs. (fO|)- 111 the general case, the endpoint 
of the critical line also depends on the RS. Instead of ( |4.5| ), which is the result for a sharp 
cutoff, we find for the general case 



M 2 



8ai 
3^2 



(6.6) 



where terms O(e) have been dropped. The entire scheme dependence is now encoded in the 
coefficient a 1; given by 



a i 



dy 



r'(y) y 



o [l + r(y)Y 



(6.7) 
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in d = 3 dimensions. This coefficient belongs to a set of expansion coefficients char- 
acterising a coarse graining scheme (see Appendix [A| for their general definition and more 
details). For each of the two classes of regulators the coefficient ai can be calculated as 



a function of the parameter n. In Fig. [L3[ the dashed line corresponds to the power-like, 
and the full line to the exponential regulator class with c = In 2. For this choice of c both 
set of regulators are normalised to r(l) = 1. For a power-like regulator, we find explicitly 
ttl = |r[l + i]r[2 - i], and for the exponential one a x = § n~ l c l / 2n {2 l / 2n - 2)T[-±]. 
It is interesting to note that although these classes of regulators do have strong qualitative 
differences, the coefficient a±, which only involves a folding of r(y) over all momenta, is 
rather stable (i.e. ±10% about the mean value). 
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FIG. 13. The expansion coefficient a\ (see text). 



We shall compare the numerical value of the tri-critical point with results obtained in the 
large -N limit via the e-expansion |SJ or a fixed dimension computation in (d = 3) B3] . As 
argued in Sect. |IIB| , the Abelian fixed point reads el = 6n 2 /N in the large -N limit, and our 
above result therefore becomes 



— 7T = loai — 

4 



N 



(9.4 



12.0) — 



(6.8) 



The e-expansion, to leading order, yields 



A, 



(54 - 136e) 



N 



(6.9) 



This is to be compared to the result of [p3j] , which reads 



A 3 96 1 1 
— = w 9 9 — 

* 7T 2 N N 



P 2 

( 3 



(6.10) 



While (|6.9| ) fails to give a reliable answer at e 
good numerical agreement with ( |6.10 ). 



1, we observe that our result ( |6.8D is in 



29 



D. Scheme dependence of the critical potential 



Here, we consider the task of computing the critical potential for coarse grainings other 
than the sharp cut-off. First, we have to obtain the corresponding fluctuation integrals. 
The most general expression (for arbitrary scheme) has been given in Appendix [B|. This 
expression still contains an integral over momenta to be performed, which is how the scheme 
dependence enters into the expression for the fluctuation integral A^. Then, the criticality 
conditions ( |4.2| ) have to be solved to find T c and p®. The sharp cut-off allowed an analytical 
computation of A k , (|B.5|) , and thus of the functions F 12 in ([Of). 



Below, in addition to the sharp cut-off, we consider the classes of power-like regulators ( |6.3| ) 
and exponential regulators (|6.4j ). From the power-like regulators, we consider the limiting 
cases n = 1 (i.e. a mass-like regulator = k 2 ) and n = oo (the sharp cutoff). As an 
intermediate case we consider also the case n = 2 (i.e. the quartic regulator Rk = k 4 /q 2 ). 
The exponential regulators are represented for n = 1 (i.e. Rk = q 2 / (exp q 2 / k 2 — 1)), and 
n = oo (the sharp cutoff). A continuity argument suggests that the critical potentials for 
intermediate values of the coarse graining parameter n should appear within those limits set 
by n = 1 , 2 and n = oo . 

No explicit analytical expressions for the coarse grained free energy have been found in these 
cases. For the mass-like and the quartic regulator we used the integrals ( |B.7| ) and ( |B.8| ), 
respectively, while ( |B.2|) is used for the exponential regulator. Then, the problem of solving 
the criticality conditions reduces to the optimization of two integral equations. 

We find that the critical temperature T c depends very weakly on the different schemes. In- 
deed, plotting T c as a function of the Higgs field mass we find that the lines corresponding 
to different schemes are almost on top of each other, inducing a relative error well below the 
1% level (and thus below the error already present due to other approximations). A similar 
situation holds for the v.e.v., where we find a relative error below a few percent. 



In Fig. [TJ|, the entire critical potential (in units of the Ad v.e.v.) is displayed for different 
coarse grainings at mn = 38 Ge V (left panel) and at mn = 70 GeV (right panel) . The labels 
s, q, m and e denote respectively the sharp cut-off, the guartic/mass-like regulator, and 
the exponential cut-off from (|6.4|) for n = 1 and c = 1. 

We first consider mn = 38 GeV, and notice that the s and q lines turn out to be on top 
of each other. Furthermore, it is realised that the v.e.v. is nearly independent on the RS, 
as is the shape of the potential close to the minima. The main dependence concerns the 
local maximum of the critical potential. This dependence will therefore affect integrated 
quantities like the surface tension, but not those related to the v.e.v., like the latent heat. 
The error for the surface tension in the present case is about a few percent. 

For mn = 70 GeV the dependence on the scheme is more pronounced than in the previous 
case. Furthermore, the v.e.v. receives - for the mass-like regulator - a sizeable shift towards 
smaller values. Again, the variance is strongest around the maximum of the critical poten- 
tial, and dominant in the non-convex region of the critical potential. The additional shift 
in the value of the v.e.v. entails a corresponding shift for the outer region of the effective 
potential, as opposed to the case for smaller Higgs field mass. 

It is interesting to make contact with the qualitative considerations presented at the be- 
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ginning of this section, and to compare the scheme dependence observed in Fig. [14] with 
the reliability of the coarse grained potential in its different regions, due to the approxi- 
mations employed. Recall that the present computation is based on neglecting the scalar 
fluctuations. This approximation is more reliable for the outer part of the potential than 
for the non-convex part of it (more precisely, around a small region of the maximum of 
the inner part of the potential). Here, scalar fluctuations cause ultimately the flattening of 
the potential in the IR limit. While we have seen in Sect. |V D| that this approximation is 
still reliable for mn = 38 GeV, we certainly expect larger corrections for mn = 70 GeV (see 
the discussion of Sects. V D and |V hj ). It is quite remarkable that the scheme dependence 
indeed seems to reflect the weakness of the approximation for this region of the potential. 
Our computation thus turns the qualitative statement into a quantitative one. 
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FIG. 14. The critical potential for me = 38 GeV (left panel) and mn = 70 GeV (right 
panel), and different regulator schemes: exponential regulator (e), sharp cut-off (s), quartic 
regulator (q) and mass- like regulator (m). 



Finally, we briefly comment on the different regulators used. It is well-known that the mass- 
like regulator is marginal in the sense that it has a poor UV behavior which makes its use for 
certain applications questionable (a more refined discussion has been given in |[22|| ). From 



Fig. [14], we learn that the critical potential as obtained for the mass-like regulator deviates 
the most from the results for the other regulators employed. Considering the class of power- 



like regulators, we see from Fig. [L4| that the width between the quartic and the sharp cut-off 
limit is significantly smaller than the deviation for the mass-like regulator. This observation 
strongly suggests that the mass-like regulator should be discarded for quantitative consider- 
ations, although it remains, in the present example, a useful regulator for studying the main 
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qualitative features of the problem.^ Discarding the mass-like regulator from our discussion, 
we end up with the observation that the error induced through the scheme dependence is 
of the same order of magnitude for algebraic as for exponential regulators. For the present 
case, and at this level of accuracy, no further qualitative differences are observed between 
the exponential regulators ( |6.4j ) and the power-like ones ( |6.3| ) for n > 2. 

In summary, we conclude that a quantitative analysis of the scheme dependence indeed 
yields non-trivial information regarding the accuracy of the approximations or truncations 
employed, as suggested by the qualitative argument presented in Sect. [VI A[ In addition, we 
have found some evidence for why a mass-like regulator, as opposed to exponential or higher 
order power-like regulators, should be discarded for accurate quantitative considerations. 
However, as the qualitative features are still well described by a mass-like regulator, and 
as the quantitative deviation is not too big, this also suggests that a mass term regulator 
could be very useful for an error estimate.^] Typically, analytical computations are largely 
simplified for mass-like regulators, allowing for a simple cross-check of the results. 



VII. SUMMARY AND OUTLOOK 

We have studied in detail the first-order phase transition of Abelian Higgs models in 3+1 
dimensions at finite temperature. Properties of the transition are determined by the un- 
derlying fixed point structure of the 3d theory such as the cross-over of the Abelian charge 
from the Gaussian to the Abelian fixed point. We computed all physical observables at 
the phase transition, the phase diagram in the domain of first-order transitions and the tri- 
critical point. The analysis has been restricted to the region of parameter space where the 
dimensional reduction scenario applies and a perturbative matching of the Ad parameters to 
the corresponding 3d ones is possible. The main contribution to the free energy (and thus 
to the physical observables at criticality) stem from the remaining effective 3d running for 
which we have used a Wilsonian renormalisation group to leading order in the derivative 
expansion, neglecting the scalar, but not the gauge field anomalous dimension. The latter 
is related to the non-trivial running of the Abelian gauge coupling, which is described by an 
effective fixed point. While this fixed point is well understood in the large -N limit where 
the tri-critical fixed point is known, its precise form is not yet established for the relevant 
case of N = 1. We therefore studied the parametric dependence of physical observables on 
the fixed point value. A quantitative discussion of the relevant physical scales, which are 
easily accessible within a Wilsonian framework, has also been given. 

The main effect on physical observables due to the presence of a non-trivial fixed point 
depends on the ratio between the cross-over scale k cr (which defines the cross-over to the 
Abelian fixed point) and the typical scales characterising the first-order phase transition 



18 This conclusion coincides with those of [22] based on more formal considerations regarding 
mass- like regulators. 

19 An error estimate based on the mass-like regulator is rather conservative as it seems to overes- 
timate the scheme dependence. 
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(like the discontinuity scale fcdi S , or Stable) ■ For k cr small as compared to Stable the observed 
dependence is weak. The sizeable deviations from the perturbative e 2 (k) ~ e 2 (A) - behaviour 
only set in at very small scales below fc sta t>ie and are no longer relevant for the phase transi- 
tion itself in this situation. The main effects are restricted to alterations in the far infra-red 
region, like the details of the flattening of the inner part of the potential. On the other hand, 
a strong dependence emerges for k cv larger than k sta .ue- 
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FIG. 15. The relative variation of T c and pQ with the effective Abelian 
fixed point. Here, the regions T/T re f and p/p re i compare the critical temper- 
ature and the v.e.v. as a function of e*. The reference values are obtained for 
e* = VGtt rj 7.7. 

Most of our results for the physical observables can be summarised as in Fig. |TJ. Here, 
the reference values T re f and p re f are given for e* = y/En (which corresponds roughly to 
k cr « fc s tabie), and for the sharp cut-off regulator. In the present approximation, the critical 
temperature is insensitive to the running gauge coupling. On the other hand, the v.e.v. ap- 
pears to be quite sensitive to the actual fixed point value, in particular for larger Higgs field 
mass. The phase transition weakens significantly for small fixed point values. The reason 
is that the gauge coupling is decreasing strongly for small fixed point values at scales larger 
than the scale where the critical potential reaches its degenerate shape, that is above the 
scale of decoupling. These results compare well with perturbation theory, except for very 
large or very small values for the Abelian fixed point. Corrections due to the non-trivial 
scaling of e 2 (k) remain below 10% for e 2 > 6tc 2 and mn below 70 GeV, but do grow large 
as soon as e 2 is a below 6n 2 . We conclude that e 2 ~ 6n 2 is a good leading order approxi- 
mation for small Higgs field mass as higher order corrections are small. For mn = 30 GeV, 
we also compared the value for the critical temperature with lattice simulations and found 
agreement below 4%. The sensitivity on e 2 < 67r 2 for larger Higgs mass, in turn, requires a 
better determination of the fixed point in this domain. This concerns in particular physical 
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observables like the critical exponents at the endpoint of the line of first order phase tran- 
sitions. 

For generic regulator function the free energy in the type-I regime has been given as an 
integral (one remaining integration). For the case of a sharp cut-off regulator, we obtained 
an explicit analytical solution for the free energy, given the non-trivial scale dependence of 
the Abelian charge. In the present article, we evaluated all relevant quantities for initial 
conditions obtained from a perturbative dimensional reduction scenario relevant for a high 
temperature (cosmological) phase transition. 

The explicit result for the effective potential can also be of use for applications to the su- 
perconducting phase transition, or for the nematic to smectic-A phase transition in certain 
liquid crystals. The main change would concern the initial potential for the effective 3d flow 
of the potential, and the numerical value of the Abelian charge at that scale. These changes 
affect in particular the ratio fc cr Astabie, an d therefore the above discussion, as both scales 
depend in a qualitatively different manner on e 2 (A) and U\. 

In addition, we studied the dependence of our results on the coarse graining procedure em- 
ployed. We have seen that the physical observables do depend only very weakly on the 
coarse graining. This is encouraging, as a strong dependence would have cast serious doubts 
on the approximations used. Furthermore, we employed a variety of qualitatively different 
coarse grainings ranging from the mass-like and other polynomial regulators over exponen- 
tial ones to the sharp cut-off regulator. Therefore, our result can be seen as an important 
consistency check of the method. The weak variation w.r.t. the coarse graining which is to 
be interpreted as an 'error bar' for the observables, is smaller or of about the same size as the 
error expected from higher order operators for the coarse grainings studied. This 'error bar' 
would vanish only if no truncation to the effective action would have to be performed. We 
also observed an intimate relationship between the truncation of the effective action, and the 
error bar introduced through the scheme dependence. More precisely, it is observed that the 
scheme dependence is largest in regions where a similarly large effect due to the neglecting 
of the scalar fluctuations in the non-convex region of the potential is expected. While this 
result is not entirely unexpected, a quantitative evidence for it has never been presented be- 
fore. It would be useful if further quantitative results in this direction could be established. 
This concerns in particular the cross-dependences between an optimal coarse-graining that 
minimizes the scheme dependence, and an optimized convergence of systematic truncations 
and approximations . 



An important open question for future work concerns the precise IR behaviour of the Abelian 
charge. This, of course, is an intrinsic problem of the 3c? theory. As argued, our current 
understanding is mainly limited due to an insufficient understanding of the field and/or 
momentum dependence of the Abelian charge. It might be fruitful to consider alternatively 



a thermal renormalisation group to improve the situation [f24]]. At the same time, the in- 
clusion of higher order corrections due to scalar fluctuations will also become important - 
close to the critical points - for a reliable determination of critical exponents and other uni- 
versal quantities. It would also be interesting to consider the S'L r (2)-Higgs theory, where a 
non-trivial endpoint of the line of first-order phase transitions has been established recently. 
A field theoretical understanding of this endpoint is still missing, and a derivation of the 
related critical indices from field theory would be desirable. Again, one expects that the IR 
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behaviour of the gauge coupling, in competition with the scalar fluctuations, is responsible 
for the existence of the endpoint. 



APPENDIX A: RS DEPENDENCE AND THRESHOLD FUNCTIONS 

The solution of the flow equation (and the related physical observables) can be written as 
momentum integrals over a measure, which depends on the precise implementation of the 
coarse graining. We employ the notation of ||13|| , where a scheme dependent measure has 
been given (in d dimensions) as 

for momentum-dependent functions f(y), where y = q 2 /k 2 , and q is the loop momenta. As 
a consequence of the conditions ( |6.1| ) on the regularisation function r(y) it follows that the 
momentum measure — r (y)/(l + r) l+d ' 2 is peaked. The measure is normalised to one, 

I r [l) = 1. (A.2) 

This implies that I r [f] depends on the coarse graining as soon as / displays a non-trivial 
dependence on momenta. 

As an example, let's consider the threshold functions in(uj), defined as 

poo r'(v)v 1+d / 2 

«M = -«*+»)/. < A - 3 ' 



They are related to the above measure through 

. 2 



pd+2 

(A.4) 



(pa + CJ )n+l 

Here, we also introduced the dimensionless effective (regularised) inverse propagator 

P 2 (y)=y + yr(y) . (A.5) 

The threshold functions can always be expanded as a Taylor series in powers of u>. Let us 
define the corresponding RS dependent expansion coefficients 

a k = I r [p- k ] , (A.6) 

which are the k th moments ofl/P w.r.t. the measure I r . These coefficients appear in the 
computation of the endpoint of the critical line (|6.6| ), which is proportional to the coefficient 
a\. For a power-like regulator r(y) = y~ n [see eq. ( |6.3|) 1 we find for arbitrary dimension d 



r[i + f + |] 

A more detailed discussion of these coefficients and a related discussion of the convergence 



a k = % r[l + fl L2 r 2 ^ (A.7) 

2 1 2ni r[i + 1 + §] y J 



of amplitude expansions and optimised coarse-graining parameters is given in [34 
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APPENDIX B: THE FLUCTUATION INTEGRAL 



The fluctuation integral reads 

1 M 

A fe (p) = 



2tt 2 



dk / dy 



2e\p r'(y) y 



5/2 



P 2 2elp + P 2 k(l + k/k c 



(B.l) 



Note that we have normalised A(0) = in the above definition. The remaining integrals in 
( |B.1|) can be solved in different ways, either first performing the momentum integration or 
the scale integration. Integrating first w.r.t. k yields (for the notation see Appendix A) 



A k (p)=I r [U(p,P)} 



(B.2) 



where 



3tt 2 U(p, P) = 2e{p / dk 



P 3 k 2 



P 2 k{\ + k/k cr ) + 2elp 



2e 2 j/P 2 + k + k 2 /k c 



= -2ei-p P *.(* - A) + 4p P ki In ^_ /pi + A + ^ 
+2elp P k cr (4elp/P 2 - k CT ) G k)A {l - 8elp/P 2 k cr ) , 
with I r defined in ( |A.1| ) and P{y) in ( |A.5| ). The function G(Q) reads 
1 fl + 2k/k CT - VVL l + 2A/k CT + Vtt s 



(B.3) 
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(B.4) 



For a sharp cut-off regulator 
analytically to give 



(fc cr + 2fc)(A; cr + 2A) 

the remaining momentum integration can be performed 



2n z A 



2 a( s ) _ 1 A 3 



±A d In 1 



2elpk c 



- ifc 3 In I 1 + 



A(A + k c 



2e 2 pk cr 
k(^k ~\~ k^ 



\kl T hi 



'A + k c 
k + k 

c: 



+ efe-6e^ cr )ln ^ (A; + fecr) + 2e ,- fec] 



+ |e^pA; cr (A - fc) 



-ifc cr (-16e^p 2 + lOe^ - k 2 cr ) G kA (l - 8e^/A; cr ) . (B.5) 
We have normalised A(p) such that A(0) = 0. 

On the other hand, performing first the scheme dependent momentum integration leaves us 
with the following remaining integrals, 
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A^(p) 
A (m) (p) 

A (9) (P) 



1 rA / 

dfc /fc 2 In 1 + 



2vr 2 7 fe 

1 M _ _ 



fc(l + k/k cr ) 



2tt A 

1 /" A r t 



2elp 



V2< 
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y\ k(l + k/k, 
Ik d * ^ (* ~ (* + 



-l/2> 



fc(l + A;/A; c 



(B.6) 
(B.7) 

(B.8) 



Here, the indices refer to the sharp (s) , the mass-like (m) and the quartic (q) cut-off function, 
as defined in Sect. |VI. 



APPENDIX C: INCLUDING SCALAR FLUCTUATIONS 



In order to obtain an estimate of the effect of the scalar fluctuations we will solve eq. (|2.9|) 
with Uk on the r.h.s. replaced by U\, for a sharp cut-off regulator. The flow equation becomes 
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,dU k (p) 



' m lApY 



kHk u V k 2 
with the masses mf given through 



k 2 
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, , 2e 2 {k)p 
k 2 



m 2 (p) 

We introduce the functions 



m 



R + A R p, m 2 (p) = m R + 3A R p . 



K{m 2 



1 f k , 9 / m 2 \ 
-4r7 2 L dyyln { 1+ f) 



127T 2 
,2, 



2m 2 k — 2m 3 arctan f — 



m 



k 6 In 1 + — 



(A; <-> A) 



J,(p) = K[mi(p)] - K[m%( 
The solution to the flow ( |C.1| ) then obtains, using also from ( B.5|) , as 

u k (p) = u A (p) + A«(p) + j x {p) + j 2 {p) . 



(C.l) 



(C.2) 



(C.3) 
(C.4) 

(C.5) 



The effect of the additional terms on the shape of the critical potential is about a few percent, 
increasing towards higher values for mn- 
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